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Recently, a strong debate has been pursued about the Newtonian limit (i.e. small velocity and 
weak field) ol fourth order gravity models. According to some authors, the Newtonian limit of f(R)- 
gravity is equivalent to the one of Brans-Dicke gravity with ujbd — 0, so that the PPN parameters 
of these models turn out to be ill defined. In this paper, we carefully discuss this point considering 
that fourth order gravity models are dynamically equivalent to the O'Hanlon Lagrangian. This is a 
special case of scalar-tensor gravity characterized only by self-interaction potential and that, in the 
Newtonian limit, this implies a non-standard behavior that cannot be compared with the usual PPN 
limit of General Relativity. The result turns out to be completely different from the one of Brans- 
Dicke theory and in particular suggests that it is misleading to consider the PPN parameters of this 
theory with ujbd = in order to characterize the homologous quantities of /(_R)-gravity. Finally 
the solutions at Newtonian level, obtained in the Jordan frame for a /(_R)-gravity, reinterpreted as 
a scalar-tensor theory, are linked to those in the Einstein frame. 
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I. INTRODUCTION 

Recently, several authors claimed that higher order theories of gravity, in particular /(i?)-gravity are charac- 
terized by an ill-defined behavior in the Newtonian regime. In a series of papers [4], it is discussed that higher order 
theories violate experimental constraints of General Relativity (GR) since a direct analogy between /(i?)-gravity and 
Brans-Dicke gravity [3j gives the Brans-Dicke characteristic parameter, in metric formalism, ujbd = while it should 
be ujbd — > 00 to recover the standard GR. Actually despite the calculation of the Newtonian limit of f(R), directly 
performed in the Jordan frame, have showed that this is not the case 0, Hj], it remains to clarify why the analogy 
with Brans-Dicke gravity seems to fail its predictions also if one is assuming f(R) ~ i? 1+e with e — > 0. The shortcom- 
ing could be overcome once the correct analogy between /(i?)-gravity and the scalar-tensor framework is taken into 
account. 

The action of the Brans-Dicke gravity, in the Jordan frame, reads : 



aBD 
A .1F 



4>R + uj B d^- 1 V X£% 



(1) 



where there is a generalized kinetic term and no potential is present. On the other hand, considering a generic function 
f(R) of the Ricci scalar R, one has : 



-9 



f(R) + XC r . 



(2) 



In both cases, X = ^5^, is the standard Newton coupling, C m is the perfect fluid matter Lagrangian and g is the 
determinant of the metric. 

As said above, /(i?)-gravity can be re-interpreted as a scalar-tensor theory by introducing a suitable scalar field 
<fi which non-minimally couples with the gravity sector. It is important to remark that such an analogy holds in a 
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formalism in which the scalar field displays no kinetic term but is characterized by means of a self-interaction potential 
which determines the dynamics (O'Hanlon Lagrangian) |6|. This consideration, therefore, implies that the scalar field 
Lagrangian, equivalent to the purely geometrical f(R) one, turns out to be different with respect to the above ordinary 
Brans-Dicke definition (TTJ). This point represents a crucial aspect of our analysis. In fact, as we will see below, such 
a difference will imply completely different results in the Newtonian limit of the two models and, consequently, the 
impossibility to compare predictions coming from the PPN approximation of Brans-Dicke models to those coming 
from /(i?)-gravity. 

The layout of paper is the following. In Sec. [Til we discuss the solutions in the Newtonian limit of /(-R)-gravity by 
using the analogies with the O'Hanlon theory. Sec lIIIl is devoted to the analysis of the solutions in the limit f(R) —> R 
and the interpretation of PPN parameters 7, j3. Conformal transformations and the solutions in the Newtonian limit 
approximation are considered in Sec lIVI Concluding remarks are drawn in Sec. [V] 

II. THE NEWTONIAN LIMIT OF /(i?)-GRAVITY BY O'HANLON THEORY 

Before starting with our analysis, let us remind that the field equations in metric formalism, coming from f(R)- 
gravity, are 



= f'(R)R flv - \f{R)g^ - f'(R)-,^ + g^nf = XT^ 



(3) 



which have to be solved together to the trace equation 



af{R)+ nR)R-mR) = * T 



Let us notice that this last expression assigns the evolution of the Ricci scalar as a dynamical quantity. Here, 

T nv = *HV g£m) j g ^ ^ energy-momentum tensor of matter, while T — T° ' is the trace, f'(R) = ^ . 

\F r 9 <V 6 ' J v ' dR 

The conventions for Ricci's tensor is R^ = R a ^av while for the Riemann tensor is R a p tlv = r2 v „ + .... 

The affine connections are the Christoffel symbols of the metric: T^g = \g^° (g aa ,i3 + gp<y,a — gap.a)- The adopted 

signature is (H ). 

On the other hand, the so-called O'Hanlon Lagrangian 0] can be written as 



<pR - V{4>) + X£ r , 



(5) 



where V(cj>) is the self-interaction potential. Field equations are obtained by varying Eq. ([5]) with respect both g^ u 
and <f> which now represent the dynamical variables. Thus, one obtains 



(6) 



R- d ^tl = 0. 



(7) 



U4>- 



, dV(cf>) 



2V{4>) 



X. 



-T , 



(8) 



where we have displayed the field equation for 4>. Eq. © is a combination of the trace of © and ([7]). /(i?)-gravity 
and O'Hanlon gravity can be mapped one into the other considering the following equivalences 



4> = f'(R) 



(9) 
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V(<f>) = f'(R)R-f{R) 



(10) 



dV{cj>) 



2V((f>) = 2f(R) - f'(R)R 



(11) 



and supposing that the Jacobian of the transformation <f) — f'(R) is non- vanishing. Henceforth we can consider, 
instead of Eqs. ©-(HI, a new set of field equations determined by the equivalence between the O'Hanlon gravity and 
the /(-R)-gravity: 



4>Rfj,u 



1 



v{4>) + <\> 



dV(cj>) 



(12) 



□</>- 



dV{cj>) 



2V( ( f>) 



-T 



(13) 



where = T^ v - ^Tg^. 

Let us, now, calculate the Newtonian limit of Eqs. (IT2|) - (IT3l) . To perform this calculation, the metric tensor 
and the scalar field 4> have to be perturbed with respect to the background. After, one has to search for solutions at 
the (v/c) 2 order in term of the metric and the scalar field entries. It is 



9^ 



(2) 

Soo 







,(2) 



(14) 



The differential operators turn out to be approximated as 



(15) 



□ sa <9n - A and 



sa d 2 



(16) 



since time derivatives increase the degree of perturbation, they can be discarded From a physical point of view, 
this position holds since Newtonian limit implies also the slow motion. 

Actually in order to simplify calculations, we can exploit the gauge freedom that is intrinsic in the metric definition. 
In particular, we can choose the harmonic gauge <7 po T^ cr = so that the components of Ricci tensor reduces to 



P( 2 ' - 1 ArP 

n 0Q — 2^00 



(2) 



R 



(3) 



(17) 



R (2) = iAq 



(2) 



Accordingly, we develop the self-interaction potential at second order. In particular, the quantities in Eqs. (|12p and 
(EGO) read : 



vw + /-m * vm + ^ dxm + L(o) fvm + 



(18) 



d 2 v((f>^) dv(4>^) 



(19) 
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Field Eqs. (1121) - (|13]) . solved at 0-th order of approximation, provide the two solutions 



Vm=0 and ^fl=0 



(20) 



which fix the 0-th order terms of the self-interaction potential; therefore we have 



(21) 



where the constant factor m 2 can be easily interpreted as a mass term as will become clearer in the following analysis 
(see also @). Now, taking into account the above simplifications, we can rewrite the field equations at the {v/c) 2 
order in the form : 



A (2) — ' - v^UJ 



2X 



i(0) 



(22) 



^ ~ 0(0)^ +m 0(0)^+^(0) 



(23) 



The scalar field solution can be easily obtained from Eq. (|24|) as : 



(24) 



^(x) =~ |d 3 x'g(x,x')T(°)(x') 



where C/(x, x') is the Green funtion of the operator A — m 2 , while, for and gff, we have 



,(2) 



(25) 



X 



2^0(0) 



^IlfflW , m 2 /- d 3 x ^ (2) (x') 



|x-x'| 47T0(°) 



(26) 



27T0(°) 



d 3 x' * J 



« ' d 3 x 



> (2) (x') 



2 



(27) 



The above three solutions are a completely general result p}. An example can make clearer the discussion. We 
can consider a fourth order gravity Lagrangian 1 of the form f(R) = aR + bR 2 so that the "dummy" scalar field 

reads = a + 2bR. The relation between <p and R is R = ^ while the self-interaction potential turns out the be 

/ / a) 2 

V(<p) = — satisfying the conditions V(a) = and V'(a) — 0. In relation to the definition of the scalar field, we 

can opportunely identify a with a constant value (f>^ — a. Furthermore, the scalar field mass can be expressed in term 



1 It is important to stress that, in the Newtonian limit of any analytic /(i?)-gravity model, we need to consider only the first two derivatives 
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of the Lagrangian parameters as m 2 — — <f>^ — '—^ — - = r- Since the Ricci scalar at lowest order (Newtonian 

3 dcp z 6b 

limit) reads 



i?~i?( 2 )(x) 



0( 2 ) (x) 



2b 



X 
66 



J dV5(x,x')T(»)(x'), 



(28) 



if we consider a point-like mass M, the energy-momentum tensor components become respectively Too = P, T ~ p 
while p = M5(x), therefore we obtain 



(27T) 1 ' 12 'r g m 2 e-™l x l 



(29) 



where r g is the Schwarzschild radius. The immediate consequence is that the solution for the scalar field </>, up to the 
second order of perturbation, is given by 



4> = a + 



(2n 



l/2 r g-m|x| 



(30) 



loo and 9if, where = 
As matter of fact, the metric solutions at the second order of perturbation are 



In the same way, one can deduce the expressions for g^' and gf^ , where Eqq = | pc 2 and £ tJ 



(0) 



kpc 2 8i 



2 ^00 °%3 ■ 



4 r„ (27T) 1 / 2 r e- m W 
5oo = l-- ' ' 



3a x 



3 a 



(31) 



9ij = - i 1 



2(2tt 



A Is 

3a |x| 

,l/2 r 



(2tt 



3 a 



- m | x | 



3 a 



- m | x | 



X{ X j 



(32) 



These quantities show that the gravitational potential coming from the O'Hanlon theory of gravity is non- Newtonian. 
The correctios have the meaning of scale parameters defining characteristic sizes and masses 0, Q • 



III. THE BEHAVIOR OF SOLUTIONS FOR f(R) -> R AND THE INTERPRETATION OF PPN 

PARAMETERS 7, /3 

The results (|31[) - (|32l) are equivalent to those obtained in /(i?)-gravity [4|. This point is very important since such 
a behaviour prevents from obtaining the standard definition of the PPN parameters as corrections to the Newtonian 
potential. As matter of fact, at the Newtonian level, it is indeed not true that a generic /(i?)-gravity model corresponds 
to a Brans-Dicke model with ujbd — 0. In particular, in such a limit, it is not correct to consider the PPN parameter 

7 = — — — (see jH) of Brans-Dicke gravity and evaluating this at ujbd = 0. In this case, one obtains 7 = 1/2 
2 + lubd „ 

as suggested in [2( and the standard Newton potential (7 = 1/2) could never be recovered. Differently, because 
of the presence of the self-interaction potential V((f>) in the O'Hanlon theory, a Yukawa-like correction appears and 
it contributes in a completely different way to the post- Newtonian limit. As matter of fact one obtains a different 
gravitational potential with respect to the ordinary Newtonian one and the fourth order corrections in term of the 
v/c ratio (PPN level) have to be evaluated in a different way. In other words, considering Brans-Dicke and O'Hanlon 
theories, despite of their similar structure, will imply completely different predictions in Newtonian limit. Such an 
achievement represents a significant argument against the claim that fourth order gravity models can be ruled out 
only on the bases of the analogy with Brans-Dicke PPN parameters. 

Another important point has to be considered. The PPN-parameters 7 and /3, in the GR context, are intended to 
parameterize the deviations from the Newtonian behaviour of the gravitational potentials. They are defined according 
to the standard Eddington metric 
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^ 2 
Ixl ' 2 Ixl 2 



<?oo = l-rfr + (33) 



r 



1+7^J%- (34) 

In particular, the PPN parameter 7 is related with the second order correction to the gravitational potential while /3 
is linked with the fourth order perturbation in v/c. Actually, if we consider the limit f(R) — > R from Eqs. pip and 
(|32|) . we have 

■900 = 1- ^-pr, (35) 
da |x| 

Since a is an arbitrary constant, in order to match the Newtonian gravitational potential of GR, we should fix a = 4/3. 
This assumption implies 

5oo = 1 - ft , (37) 

*«=-( i+ 5rK (38) 

which suggest that the PPN parameter 7, in this limit, is 1/2, that is in striking contrast with GR predictions. 
Such a result is not surprising. In fact, the GR limit of the O'Hanlon theory requires <fi ~ const V(4>) — > but 
such approximations induce mathematical inconsistencies in the field equations of /(i?)-gravity, once these have been 
obtained by a general O'Hanlon Lagrangian. Actually, this is a general feature of the O'Hanlon theory. In fact it can 
be easily demonstrated that the field Eqs. (fT21) and (fTU|) do not reduce to the standard GR ones since we have : 

X 

R^iv — j (39) 

jT = (40) 

but Yu^ v components read £00 = \p and £y = \p&ij = |£oo<% instead of Sqq = \p and S%j = \p5ij = SooSij as 
usual. Of course, = T M „ — \g^ v T with T^ u the energy-momentum tensor of matter and Einstein equations are 
written in the form 



R^v — ^ Sfj, u . (41) 

Such a pathology is in order even when the GR limit is recovered from the Brans-Dicke theory. In such a case, in 
order to match the Hilbcrt-Einstein Lagrangian, one needs <p ~ const and lj — 0, the immediate consequence is that 
the PPN parameter 7 turns out to be 1/2, while it is well known that Brans-Dicke model fulfils the low energy limit 
prescriptions of GR in the limit uj — > 00. Even in this case, the problem, with respect to the GR prediction, is that 
the GR limit of the model introduces inconsistencies in the field equations. In other words, it is not possible to impose 
the same transformation which leads BD theory into GR at the Lagrangian level on the solutions and the observables 
obtained by solving field equations descending from the general Lagrangian. The relevant aspect of this discussion 
is that considering a /(i?)-model, in analogy with the O'Hanlon theory and then supposing that the self-interaction 
potential is negligible, introduces a pathological behaviour on the solutions leading to the PPN parameter 7 = 1/2. 
This is what happens when an effective approximation scheme is introduced in the field equations in order to calculate 
the weak field limit of fourth order gravity by means of Brans-Dicke model. Such a result seems, from another point 
of view, to enforce the claim that fourth order gravity models have to be carefully investigated in this limit and their 
analogy with scalar-tensor gravity should be considered accordingly. 
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IV. THE CONFORMAL TRANSFORMATIONS AD THE NEWTONIAN LIMIT 



Along the paper we have discussed the Newtonian limit of f{R)- gravity in term of scalar-tensor gravity rigorously 
remaining in the Jordan frame. In this section, we discuss the Newtonian limit when a conformal transformation 
is applied to the O'Hanlon theory. In other words we discuss fourth order gravity models in the Einstein frame in 
place of the Jordan one when a redefinition of the metric in the conformal sense is performed. A scalar-tensor gravity 
theory is, in some sense, a generalization of both the Brans-Dicke and the O'Hanlon theories, that is 



A ST 



F(<f>)R 



V{4>) + X£ r . 



(42) 



Such a theory can be transformed by means of a conformal transformation g^ v = A(x)g M „, with A(x) > satisfying 
the condition F(^)A(x) _1 = A e E, as 



rf 4 - 



A% F = I "'•■"V ".'/ 
The relations between the quantities in the two frames are 

n(ip)dip 2 = a 



AR + n(ip)ip, a ip'' a - W{lp) + xc r , 



(43) 



"(4>) _ 3 / dlnF (4>) 
F(4>) 2 



d(j) 2 



(44) 



r ^r, 



F(4>(.¥>)) 



In the case of the O'Hanlon Lagrangian, (O, i.e. F((f>) — <fi, w(0) = 0, the previous Lagrangian turns out to be 
simplified and the transformation rule between the two scalar fields reads 



0(p)V = — ^ => ^fce ± ^V /7 W^, (45) 
2 q> £ 

where k is a integration constant. If we suppose fl(ip) — — ilo < 0, we have 



ke ±Y *. 



(46) 



where Y = 



The transformed action ([5]) in the Einstein frame is 



aOH 
■^EF 



efW- 



A 2 y»2 

AR- n <P;af ;a - ^e^V(ke ±Y n + ^-e^C r , 



A 



(47) 



The field equations are now 



' 2n n<p - £ e T2Y V[ §v [ke ±Y V) -p 2YV(ke ±Y ^)} + X^ = 
where the matter tensor, which now coupled with the scalar field <p, in the Einstein frame |l4j reads 
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rpip 



-2 5{^gC m ) 

k 2 ' 



T 2Y V 



A 



e TYip g P . 



(49) 



and 



dip 



A 2 Y 



k 2 



If 



8 £ m / A 



e TYip g P c 



(50) 



Actually, in order to calculate the Newtonian limit of the model in the Einstein frame, we can develop the two scalar 
fields at the second order <j> ~ <j)(°> + 4>^ and <p ~ tp^ + ip( 2 > with respect to a background value. This choice gives 
the relation : 



(0) = ± y-l ln : 



w (2) _ ±y-iCr 



Let us consider the conformal transformation 
if </>(°) = A, that 



(51) 



ig^v From this relation and considering the (|46p one obtains, 



( -(2) _ (2) </,< 2 > 

5oo — 5oo "+" 



(52) 



Sllx. . 

0(0) 



As matter of fact, since g$ = 2$ JF , g\f = 2* JF 5y and g$ = 2<P EF , gjf = 2^ EF <5 y from {52J| it descends a 
relevant relation which links the gravitational potentials of Jordan and Einstein frame : 



^EF 



0< 2 > 
20C) 

20(°) 



(53) 



If we introduce the variations of two potentials: A$ 
relevant result of this section: 



§JF _ $EF and ^ = qJF _ ■q/EF we b tam tne most 



Ad> = -A* = = T-^ 2 > « & ex /"(# = 0). (54) 

From the above expressions, one can notice that there is an evident difference between the behavior of the two 
gravitational potentials in the two frames. Such achievement suggests that, at the Newtonian level, it is possible to 
discriminate between the two frames. Specifically, once the gravitational potential is calculated in the Jordan frame 
and the dynamical evolution of cj> is taken into account at the suitable perturbation level, these can be substituted in 
the first of Eqs. (|53l) to obtain its evolution in the Einstein frame. The final step is that the two potentials have to 
be matched with experimental data in order to investigate what the physical solution is. 



V. CONCLUSIONS 



In this paper, we have used the analogy between the /(i?)-gravity and the O'Hanlon theory to discuss, in the Jordan 
frame, the Newtonian limit of the theory. The main result is that it is not possible to consider the analogy between 
/(-R)-gravity and the Brans-Dicke theory to achieve the correct PPN limit, as done several times in literature, since 
the result lobd = implying 7 = 1/2 is a pathology of the theory (both f{R) and Brans-Dicke). This means that 
the PPN-parameters have to be redefined accordingly in the Jordan frame of /(i?)-gravity without transforming the 
theory. 
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When we perform the GR limit of /(i?)-gravity, the correspondence is needed at any level (Lagrangian, field 
equations and solutions) between the GR and the /(-R)-gravity, in particular as soon as f{R) ~ R 1+e with e — » 0. 
This means that the statement 7 = 1/2 for any f(R) [2] is not correct since the transformation in terms of Brans- 
Dicke theory does not work . Supposing to modify the Hilbert-Einstein Lagrangian, the correction to the solutions 
can not produce the same displacement from Schwarzschild solution for any /(i?)-gravity and such displacement could 
not be independent from the analytical form of f(R). In fact when f(R) — > R the solutions of field equations are 
not the solutions of GR. Furthermore the Eddington parameterization [T(| is based on the hypothesis that metric has 
to match second order differential field equations. This means that gravitational potential admits the same Green 
function of the Newtonian theory [HI). In /(i?.)-gravity case, field equations are fourth-order in metric approach. 
The field equations, in Newtonian limit, admit Yukawa-like corrections and the standard Eddington parameterization 
cannot work [il. [T2I [l3j. On the other hand, in order to compare results in Einstein frame and Jordan frame, one can 
perform the Newtonian limit in both frames and then compare the solutions. Immediately it emerges that results are 
different. As final remark, it is worth saying that /(i?)-gravity can fully evade Solar system tests but results have to 
be carefully analyzed in the right frame. Forthcoming experiments could clearly give indications in this sense [TH ]. 
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